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ABSTRACT 
The Navy’s Impact Burial Prediction Model creates a two-dimensional time 


history of a bottom mine as it falls through air, water, and sediment. The output of the 
model is the predicted burial depth of the mine in the sediment in meters, as well as 
height, area and volume protruding. Model input consists of environmental parameters 
and mine characteristics, as well as parameters describing the mine’s release. The model 
user seldom knows many of these parameters, and those that are known may be of 
questionable precision. 

In order to determine which parameters had the greatest effect on the model and 
which could be simplified or eliminated, a series of sensitivity tests were performed. It 
was found that the model data ingest could be greatly simplified without sacrificing 
accuracy too much. However, several parameters including sediment shear strength were 


found to have a large effect on the model and were investigated further. 
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I. INTRODUCTION 


In minesweeping, success often hinges on knowing-as much as possible about the 
mines that have been placed and the effects the environment has had on that placement. 
Since bottom mines cannot be searched for visually, and are often difficult to locate with 
conventional sonar, an estimate of the area or height of the mine protruding from the 
sediment, or the burial depth if completely covered, is crucial information for the 
planning and execution of mine clearance operations. Determining the likely mine burial 
depth requires numerical models of the burial process and knowledge of the environment, 
including sediment properties, waves, tides, and water depth. 

Ocean deployed mines currently used by the U.S. and other nations fall into three 
general categories: bottom mines, moored mines and drifting mines. Bottom mines rest 
on the ocean floor and are generally deployed in littoral regions. Common placements | 
for bottom mines include shipping channels, harbors, anchorages, rivers and estuaries. 
Bottom mines are deployed in one of three ways: aircraft, surface ship or submarine. 
Although mines are designed to be deployed by a specific platform, most mines can be 
deployed by surface ship with little modification (NMWEA, 1991). 

Several numerical models have been developed to simulate the mine burial 
process, and constitute the only viable means for determining a predicted burial depth, 
which is critical information when clearing an area of mines. The Impact Burial (IB) 
model was developed to determine the depth at which the mine comes to rest in the 
sediment upon impact. Originally created by Arnone and Bowen (1980) at the Naval 


Coastal Systems Center in Panama City, Florida, the IB model was designed to create a 





two-dimensional time history of a cylindrical mine as it falls through air, water, and 


sediment phases (Fig. 1). The burial depth of the mine in the marine sediment is then 


é 


calculated from the mine’s velocity on contact with the sediment and the sediment 


characteristics. 


Several revisions and changes have been made to the model to refine the physics 


and allow for more realistic geometry and more extensive input from the user. Most 


notable are the changes made by Satkowiak (1987) and Hurst (1991). Other revisions 


implifying the data entry process. 
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involved translating to newer computer languages and s 
‘Water Depth and 





sediment. Labels on the right are parameters used by the model to calculate velocity, attitude, 


Figure 1. The trajectory of a cylindrical mine as it falls through three phases: air, water and 
and burial depth of the cylinder. 





Currently, the model allows the user to input nearly any value for each 
environmental parameter. Many of these parameters are rarely if ever known by the 
technician, and their inclusion makes it very difficult for the field user to get an accurate 
solution. With this in mind, a sensitivity test was designed and executed with the 
objective of simplifying the input parameters without compromising thé accuracy of the 
model’s output. By determining which can be eliminated or simplified without 
sacrificing accuracy, a complex model may be made more manageable for the customer. . 
Several parameters were found to lend themselves well to simplification, while others 
had so little effect on the outcome that the default value is sufficient. Some parameters, 
however, have such a large influence on the final burial depth and are known with such 
poor precision that a method for more accurate field measurements should be 


investigated. 








Ii. DYNAMICAL ANALYSIS OF A 
CYLINDRICAL MINE 


A. GEOMETRY OF A CYLINDRICAL MINE 
As a cylinder falls through a fluid medium it has a varying angle with respect to 


the vertical called its attitude, ® (Fig. 2). The attitude of the cylinder determines the area 
of the cylinder that is perpendicular to the horizontal, called the projected area, Ay. As 
the attitude varies with time, the projected area also changes, as do the magnitudes of the 


vertical and horizontal forces acting on the cylinder. 





Figure 2. Horizontal length, L,, as a function of cylinder diameter and length, D 
and L, and attitude ®. Rotation rate, m, determines the attitude at a given time. 





B. ATTITUDE OF A CYLINDRICAL MINE 
All forces and cavity calculations are affected by the attitude of the cylinder, ®. 


Forces are broken down into those that act parallel to the axis of the cylinder, axial 
forces, and those that are perpendicular to the axis, called cross forces. If required, these 
force values can then be geometrically translated into horizontal and vertical forces. 

If the mass is uniformly distributed within a mine the center of gravity, cg, and 
the center of buoyancy, cb, are the same. If they differ, the moments about the center of 
buoyancy which result from forces in the horizontal and vertical directions are denoted as 
Xegcb aNd Zeg-cp. These values are used along with the weight, W, and velocity, V, of the 
mine and the kinematic viscosity and density of the fluid, A, and p to compute steady 
state attitude ©" for a specific time increment (Amone and Bowen, 1980) 


or-| -\A, pV (Cosby at 


-\A; pV (Sinoy ee (1) 


C. ADDED MASS 


Following a mine's deployment by aircraft, it penetrates the air, water, and 
sediment. The periods during which the mine falls through each of these media are called 
the air phase, water phase, and sediment phase. A solid body moving in a fluid originally 
at rest behaves like a body of increased inertia (Mises, 1959), depicted using the term 
added mass, Maddea. The summation of the added mass and the mass in air, M,, is called 
the resultant mass, M;: 

M = Mz + Maddea (2) 

Where M, is the weight in air (M,=W/g), W is the weight of the cylinder and g is the 


gravitational acceleration. 








The added mass, also called apparent mass, is a function of the shape of the 
object. It is computed using inertial coefficients that account for the flow processes about 
the object and the cylinder volume, B. For a cylindrical object two inertial coefficients 
are required, K, and K2, to account for axial flow and cross flow: | 

M, = Kk} p BCos ®+ kK, p B Sin ® (3) 

D. PROJECTED AREA 
The projected area of a cylindrical mine varies as a function of the cylinder's 


attitude (Fig. 3). The horizontal and vertical lengths of a cylindrical mine, L, and Ly, are 


computed by: 
Li = L Sin® + D Cos® 
| a 


L Cos® + D Sin® (4) 


The projected area for a horizontal cylindrical mine, A, where ® = 90°, or for a 
vertical mine, A, where ® = 0°, is defined as: 


An mr LD 
A, = BD7/4 


Thus, the projected area for a cylindrical mine with an attitude of ® is computed by: 


A = Ay Sin* + A, Cos” ® (5) 











Figure 3. Parameters used to compute projected lengths and projected area. 


E. DRAG COEFFICIENT 


Drag on a circular cylinder is proportional to the sum of forces acting both 
vertically and horizontally. The drag coefficient used to calculate the total drag must 
therefore take each force into account. The drag coefficient, Cp, is derived from the 
Reynolds number and the shape and dimensions of the object. For a cylinder, the method 
used depends on whether it is a right circular cylinder or if it has a taper or rounded nose. 
Both cases are described below. 

The drag force on a cylinder is computed using the projected area of the immersed 
body, A;, the velocity of the object, V, and the drag coefficient (Schlichting, 1979): 


FK, = 0O5CpAV'A (6) 











1. Drag Coefficient for a Circular Cylindrical Mine 
In order to determine Cp we must first determine the Reynolds number for the 


object. A smooth right circular cylinder has a Reynolds number of less than 10°, 
indicating that laminar flow is dominant. The Reynolds number is calculated separately 
for axial and cross flow using the corresponding velocity, V, and the kinematic viscosity 
of the fluid, Ay: | | 

Re = VD/A, (7) 

The drag coefficient for the skin friction, Cr, is computed using the Reynolds 


number (Mises, 1959), which is then used to calculate the base drag coefficient, C*p: 


Ce 1.328 (Re) (8) 


C*p 0.33D/L + Ce [(3L/D) + 3(D/L)*} (9) 

Further modification accounts for the projected surface area, which will always be 
less than or equal to the projected area of a horizontal cylinder. This is done by applying 
a ratio of A, and the standard wetted area of a cylinder (equal to projected area of a 
horizontal cylinder), Ay. A generic coefficient of 1.1 is also applied to account for 


surface imperfections on the cylinder (Arnone and Bowen, 1980): 


Cp =  Cp*11(As/Ay) (10) 


2. Drag Coefficient for a Tapered or Round-Nosed Cylindrical Mine 
If the circular cylinder has a rounded nose or is tapered slightly, turbulent flow 


will be present and the Reynolds number is computed by the Prandel-Schlichting 


(Schlichting, 1979) skin-friction coefficient: 


Crp =  0.455/(log Re)*®* (11) 








The overall drag coefficient is calculated using both the diameter of the base, Dp, and the 
maximum diameter, D (Arnone and Bowen, 1980): 


Cp’ 
Cp 


Cr (Aw/As) * [1 + (60(L/D)? *) + 0.0025(L/D)} 
Cp’ + (0.029 (Dp/D)’ (Cp’)*”) (12) 
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Ill. MOMENTUM BALANCE OF A CYLINDRICAL 
MINE IN AIR AND WATER 


A. EQUATION OF MOTION 


The movement of a mine in the air and the water is depicted by the momentum 
equation: 

M, dV/dt = Fy.t+ F,+ Fy (13) 

Where V is the velocity of the mine, F,,,, is the force due to the air weight of the mine, F, 
is buoyancy force and Fyis drag force from equation (6). 

Buoyancy force is the upward force exerted upon a mine in the gravitational field 
by virtue of the density difference between the mine and that of the surrounding fluid. We 
use the Cartesian coordinate system (x, y, z) with the z-axis in the vertical direction, and. 
use the unit vector k along the z-axis (pointing downward). The buoyancy force is then 
computed using the density value for air or water, p, and the volume of the mine, B: 

Fh= -pgBk (14) 

B. MINE ENTRY INTO WATER | 

A cylindrical mine penetrating into water passes through two distinct regimes. 
The first regime is the cavity regime. As the mine pushes into the air-water interface, it 
creates a cavity that consists of a combination of air and water particles. The ratio of air 
to water in the cavity decreases until the fluid properties become that of water only, at 
which time the mine is in the fully wetted regime. A temporal variation of the mine’s 
vertical position can be calculated using a method similar to that depicted in the previous 


section. 


il 





1. Cavity Regime 
Upon impacting the water, the cylindrical mine enters the cavity regime of the 


water phase. The fluid in the cavity is actually a combination of air and water, and the 
properties change continuously with the mine's downward vertical motion. 

The most critical properties of a fluid when predicting the movement of an object 
within the fluid are density and kinematic viscosity. In the cavity regime, the fluid 
density is a gradient from that of air, 1.29 kg/m’, to that of water, 1025 kg/m>. In order to 
describe the ratio of water volume to total cavity volume we use the term void ratio, fp. 
This value is used to determine the actual density of the fluid in the cavity. Although it is 
difficult to determine fp precisely, the gradient is assumed to be linear. When the mine is 
fully in the air the void ratio is zero, and when fully in water it is one. The void ratio is 
modeled based on the assumption that the ratio of water volume to air volume in the 
cavity increases with an increase in the ratio of hydrostatic pressure to dynamic pressure. 
When the pressure generated by the dynamic loading of the cylindrical mine as it impacts 
the water, called dynamic pressure, is balanced by the hydrostatic pressure of the water, 
Ps, the cavity collapses and the mine 1s considered fully wetted (Arnone and Bowen, 
1980). The dynamic pressure, Pp, is the summation of the atmospheric pressure at the 
ocean surface (assumed constant) and the pressure caused by the mine movement: 

Pp = P. + % Area V" (15) 

Where Arry is the two-phase mixture density in the cavity. 
Arnone and Bowen (1980) further demonstrate this relationship by stating that the 


void ratio equals the square root of the pressure ratio: 


to = Jf F,/Ps a) 


12 





The two-phase mixture density, Array, and kinematic viscosity, Arg, in the cavity 
are computed at each time step. The fluid properties of density, p, and py. and kinematic 
viscosity, A,and Ay, for air and water are also required: 


Arm = A, (1 — fo) + foAw | (17) 
ATEM [(As Ag) (Aw Aw)] [1 — fo) Aa Aa + fo Aw Aw] (18) 


At the each time step, the hydrostatic and dynamic pressures are calculated from 
the velocity of the cylindrical mine, and the void ratio is determined. When the 
hydrostatic and dynamic pressures are equal, the void ratio takes the vale of 1 and the 
cavity collapses. The mine is then considered to be in the fully wetted regime. At this 
point the velocity of the mine is called the exit cavity velocity. If the actual water depth 
is less than the depth at which the cavity collapses, the mine impacts the sediment prior to 


exiting the cavity and there is no fully wetted regime. 


2. Fully Wetted Regime 
The trajectory of a mine in the fully wetted regime is similar to that in the air 


except the properties of the fluid have changed. The exit cavity velocity is used as an 
initial condition for the mine momentum equation (13), used to determine the mine 
trajectory. When the vertical distance of the mine traveling in the water equals the water 
. depth, the mine velocity is called the bottom impact velocity, which is the initial 


condition for determining the mine burial depth in the sediment. 


C. MINE BURIAL IN SEDIMENT 


1. Viscosity and Shear Strength 
Penetration of the cylindrical mine into the bottom sediment depends primarily on 


the attitude and velocity of the mine upon impact, as well as the sediment properties of 


density and shear strength. Initial impact of the cylindrical mine into the sediment 
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creates a cavity in which the fluid properties of water and sediment are interacting. To 
determine the mine burial depth in the sediment, one must determine the void ratio and 
the two-phase density and kinematic viscosity. A procedure similar to that depicted in 
the previous section can be used to determine these parameters. However, the 
computation of the two-phase kinematic viscosity for the water-sediment cavity differs a 
bit from that of the air-water cavity. The kinematic viscosity of the sediment, As, is not a 
pure constant, but rather is equal to the water viscosity, Ay, plus that resulting from the 
shear stress of the sediment: 
Ag = Aw + Su/ (A; dV/dz) (19) 


Where A, is the density of the sediment and S, is the shear strength. 


2. Mine Burial Dynamics 
The vertical momentum balance of a mine in the sediment phase is given by: 


M, dV/dt = Fy. + Fy + Fa+ F.+ F, (20) 
where F,, is the buoyancy force in the sediment, F, is the somes force, and Fs is the 
shear force. F, and F, are additional forces (different from air and water phases) exerted 
on the mine by the sediment. They are proportional to shear strength of the sediment and 
the projected area of the mine. If the mine is a right circular cylinder, the compressive 
force is twice the shear force: 


F.=2F. 
F,=S,A (21) 


The mine burial depth is predicted by integrating (20) with respect to time until 
the mine velocity becomes zero. Accurate values for sediment properties are essential to 
the accuracy of this process. Shear strength and density have a strong impact on the 


computation of all forces as well as buoyancy weight and added mass. 
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IV. MINE IMPACT BURIAL PREDICTION MODEL 


A. DEVELOPMENT HISTORY 
Based on the fluid dynamics of a cylindrical mine, Arnone and Bowen (1980) 


developed a mine impact burial prediction model called the Impact Burial Prediction (IB) 
Model. The basic model created a two-dimensional free-fall history of a right circular 
cylinder as is fell through air, water and sediment phases and calculated the burial depth 
of the mine after it came to rest in the sediment. The initial model had a number of 
shortcomings, particularly when handling soft or hard sediments or certain environmental 
conditions. It also handled the water-sediment cavity regime poorly. 

Satkowiak (1987) modified the IB model, including: 
Corrected reference flow used in drag calculations 
Corrected added mass term equation 
Reworked equations for sediment-cavity regime 
Allowing for drag due to cylinder’s nose 
Allowing for rounded noses 


Inclusion of water temperature 
Inclusion of retarding forces in semi-solid sediment 


The IB model was extensively revised by Hurst (1991). Equations for forces 
acting on the mine were redefined and equations to simulate rotational movement of the 
mine were added. Five primary areas were addressed and refined: 


Calculation of fluid drag 

Calculation of forces at air-water interface 

Calculation of forces upon impact with the sediment 

Allowing for rotational movement of mine as it falls 

Allowing for multilayered sediments (formerly only the deepest point was 
considered) | 


Also added in 1991 was the ability of the model to calculate fall angle dynamically, as it 


would in reality. 
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In 1993, Mulhearn’s formulation for sediment bearing strength (1993), which 
takes sediment shear strength as well as the object’s mass and geometry into account, was 
added to the IB model. The model was also modified at this time to allow for displaced 
centers of mass. 

Currently, the twenty-sixth edition of the IB model is available in ANSI C format. 
The model features interactive data input with the user being prompted at every step. No 


user manual is required, and none exists at this time. 


B. IB MODEL INPUT PARAMETERS 


The IB model has a variety of parameters that must be entered by the user prior to 
running the model. Several of these parameters are unknowns, even to the mine 
sweepers, and others are easily determined based on intelligence and published charts or 
atlases. 

One possible unknown parameter that is required is mine type. If the specific 
type of mine is known, the user enters the mine’s characteristics of air weight, water 
weight, length, maximum diameter, taper length, base diameter, and center of mass offset 
distance from the published mine lists. If unknown, an educated guess must be made and 
a profile selected from those available in the model (Table 1). Most of the mine profiles 
available in the model, and most mines used today, are right cylinders with a center of 
mass in the geometric center of the mine. The model also has the capability to handle 
tapered mines and offset centers of mass. Both the weight and geometry of the mine have 


a large effect on the model’s output, so care should be taken in selecting the mine type 


(Fig. 4). 
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Table 1: Parameters for mines that are available for selection in the model. Center of mass offset 
difference can also be specified (in meters) but is 0 for all selections available in the IB model. 
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The initial orientation of the mine can be set between 0° and 90°, with 0° being a 
vertical initial orientation and 90° being horizontal. This parameter may be based on the 
knowledge of how the mines were laid, but if unknown the default of 90° will yield the 
minimum burial depth. Another unknown value is the rotation rate, the constant rate a 
mine would rotate in degrees per second. Because rotation rate is known only in 
theoretical cases, the default has been set at zero. 

If dropped vertically with no initial rotation rate, the mine will still tend to acquire 
a horizontal velocity component as it falls through the air and water phases. This results 
in the mine impacting the sediment at a random angle, significantly altering the actual 
burial depth. Satkowiak (1988) suggested that since the attitude of the mine upon impact 
with the sediment cannot possibly be known, even for controlled tests, the model should 
be used mainly to determine the maximum and minimum limits of the burial depth. Since 
0° initial orientation provides the shallowest burial and 0° the deepest, the limits can be 


easily determined and the range of burial depth used with confidence. 


C. AIR AND WATER PHASES 


As the mine passes through the air and water phases, all forces acting upon it are 
calculated and summed. The acceleration is then calculated and integrated to provide the 
velocity of the mine as it enters into the next phase. The calculations for the air and 
water phases are similar, and for the most part are handled by a generic fluid subroutine. 
The first thing computed at each time step by the model is the effective axial and cross 
masses in air, cross and axtal velocity, and the distance from the bottom of the mine to 
the center of mass and center of gravity for the current attitude of the mine. These values 


are then used to calculate the drag coefficients, drag force and torque on the mine. 
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Finally, the velocities are translated to x and z movement (in the Cartesian 
coordinate system) and ® is recomputed to represent ®*, the steady state attitude at the 
end of the time increment. All values are calculated for a time increment of 0.01 
seconds, and the values are printed to the screen in increments of 0.5 seconds. 

The interface between each phase has a cavity regime where the nose of the mine is in a 
cavity that it actually a mixture of the two phases. The composition of the mixture _ 
changes rapidly as the mine pushes through the cavity. Depth of the cavity is calculated 
from the previous time step’s values for attitude and depth of mine. The cavity exists 
until dynamic pressure equals hydrostatic pressure, at which point the cavity collapses. 
D. SEDIMENT PHASE 

The physical properties of the sediment used by the model are the depth, density 
and shear strength of each layer. This information may be available for some areas, but 
no database currently exists that contains sediment profiles including shear strength. This 
value is independent of density in saturated sediment and 1s difficult to measure, 
especially for soft sediments where a sample profile is extremely difficult to obtain. 
There are three sediment profiles available in the model that represent hard, medium and 
soft sediment types. All mines and all initial conditions produce complete burial in the 
soft sediment case. The most variance in burial depth comes from the medium-density 
sediment profile, and for this reason it was used in this sensitivity study unless otherwise 
noted. 

The main contributions on the slowing of the _— as it impacts the sediment, 
according to Hurst, are bearing strength of the sediment (70%), hydrodynamic drag of the 


sediment (25%), and buoyancy in the sediment (5%). 
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1. Bearing Strength 
In the original IB model, bearing strength was assumed to be the undrained shear 


strength times ten. It was soon determined that this simple assumption caused the burial 
depth to be oe underestimated. By using a method devised by Mulhearn (1993) 
for calculating the bearing strength that incorporates the mine dimensions and the impact 
velocity, a more accurate result is obtained. Mulhearn bearing strength, BSy, is 
calculated using sediment shear strength, S,, and the rate of loading factor, rif: 

if = (0.5 (mine impact velocity) / D)° 

BSmu = 5.14 Sy (1 + D/ (5.1414) (1 + 0.4z/D) (rf) (22) 
Where Ly is the horizontal projection of mine, D is the mine’s diameter and z the depth 
from the mine’s lower surface to the water-sediment interface. 

In the IB model, this formulation can be toggled on or off by the user. Best 
agreement between measured and modeled burial depth was found when this formulation 
was used and shear strength was assumed to be constant for depths greater than 0.55 


meters. The correct formulation for the rate of loading factor is still uncertain, however. * 


2. Hydrodynamic Drag in Sediment 
Hydrodynamic drag in the sediment is the force required to push the sediment 


aside. Assuming a low Reynold’s number of 10°, accurate for a right circular cylinder, 
and using a standard cavitation correction of 0.55, a drag coefficient of 0.98 can be 
empirically derived. Cpy, the hydrodynamic drag coefficient, is then calculated for the 
vertical velocity, C,, at each time increment: 

Coy = (0.98 * C,7)/2 (23) 
Cpx and the calculated area in contact with the sediment are then used to calculate the 


force this drag has on the cylinder. 
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3. Buoyancy in Sediment 
Buoyancy in the sediment is due to the cavity formed as the nose pushes against 


the sediment and forms a cavity. Buoyancy is a result of hydrostatic pressure and is 
calculated using the density of the sediment p,, density of the mine p,,, and the depth in 
the sediment h: 

Fo = g h (Ps - Pm) | (24) 
It is assumed that the hydrostatic forces work uniformly on the surface area of the leading 
edge of the mine that is in contact with the sediment. The total buoyancy force is then 


calculated using the surface area of the leading edge and the buoyancy coefficient above. 


E. ENVIRONMENTAL PARAMETERS IN THE IB MODEL 


The altitude from which the mine is released determines the velocity and attitude 
of the mine as it reaches the air-water interface. If a mine does not fall straight down but 
rather “tumbles” with a constant rate of rotation, simulated in the model by providing a 
rotation rate 6, the attitude of the mine upon reaching the water is impacted greatly by the 
release altitude. Although not accounted for in the model, this rotation rate may be 
caused or affected by wind. 

In the water phase, this rotation rate is damped significantly. However, it still has 
a great effect on the angle the mine makes with the sediment upon impact. Currents may 
affect the rotation rate in the model, but again are not accounted for in the model. The 
water depth only has an effect on impact velocity if it is less than that required for the 
mine to reach terminal velocity, the velocity at which the deceleration due to frictional 
_ drag is equal to the acceleration from gravity. The velocity at aa this a is 
reached is a function of the weight of the mine. Since mines are laid in shipping channels 


almost exclusively, one may assume that water depths in excess of that required for a 
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mine to reach terminal velocity are the norm. Water temperature has an effect a the 
viscosity of seawater, and hence increases the drag of the seawater on the mine. 
Properties of the sediment are represented by density and shear strength profiles. 
Density of marine sediment tends to have a s-shaped profile with sharper gradients as 
density increases. Shear strength, the ability of the sediment to withstand pressure 
without deforming, also typically has a s-shaped profile and increases with distance from 
the water interface. The shear strength is related to the level of cohesion between the 
sediment particles. The density range of concern to the mine impact burial problem is 
1375 to 1600 kg/m’. Factors contributing to shear strength are the type of material, water 
content, history of stress or disturbance and time since deposition (Noorany, 1985). 
Although both increase with distance from the interface, there is no clear correlation 
between shear strength and density. Figure 5 is a scatter plot of density and shear 
strength values for 62 sediment samples, all taken at the water-sediment interface. For 
this particular data set, the correlation is extremely weak. Shear strength at the water- 
sediment interface can be measured in situ with a vane penetrometer or other instrument. 
A profile of shear strength such as is called for in the model must be measured 
from a core sample in the laboratory. This process is time consuming and expensive, and 
no database of shear strength values currently exists. The term bearing strength, as used 
in the IB model data ingest, refers to the undrained shear strength times 10. This value, 
however, 1s converted back to shear strength and used in the Mulhearn bearing strength 


equation. 
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Although not a parameter considered in the IB model, wave action has a direct 
effect on water depth and, therefore, on velocity of the mine as it reaches the sediment 
interface. This effect only becomes significant when the ratio of water depth to wave 


height is high, and only at very low release altitudes (Table 2). 
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Table 2. Effect of wave height on burial depth. Wave heights were set around a mean water 
depth and release altitude was selected to represent a ship or aircraft delivery. If released from 
150 meters, wave height has little effect. If released from near the water surface, wave height 
upon release has a significant impact if the depth/wave height ration is high. 
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Figure 4. Sensitivity of IB predicted burial depth (m) to various mine 
types. Characteristics of mine types are in table 1. 
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V. SENSITIVITY STUDY FOR THE 
CYLINDRICAL MINE 


The purpose of this study is to ascertain which variables the model is most 
sensitive to and which can be simplified or eliminated in order to simplify its use. Since 
some variables are typically unknown by the user, it is important to determine which of 
| these have the most impact on the model and which can be reduced to toggled values or 
default values without greatly eosin the accuracy of the model. The model was 
altered to allow most parameters to be set and a loop run of one variable at a fe The 
range of each variable was set to represent all possible conditions the model would be 
used = It should be noted that wind and currents are not accounted for in this model. 
However, the only impact they would have would be on the attitude of the mine as it 
enters each phase. All runs were made with preset mine profile “Korean Mine”, which 
has a dry weight of 538 kg, a wet weight of 251 kg, and a uniform diameter of .475 
meters. 

Since the model calculates burial depth and then geometrically calculates the height, 
volume and area protruding, these values are proportional. To confirm this, we created 
derivative plots of these values for one case and found the shapes of the curves to be very 
similar (Fig. 6). Burial depth is used to explain most of the sensitivity test, except where 
height protruding is more descriptive. 

A. SENSITIVITY TO RELEASE MEDIUM PARAMETERS 

Figure 7 demonstrates the variation of the release medium parameters of altitude, 
water depth and water temperature. Altitude, when varied from 0 to 1000 meters, has a 
small impact on burial depth (relative difference of 18%). When a more realistic upper 


limit of 300 meters for a mine laying aircraft is applied, the relative difference drops to 








just 9%. Water depth has an effect on the burial depth only if less than the depth needed 
for a mine to reach terminal velocity, in this case 20 meters. At depths greater than this 
value, the mine reaches terminal velocity in the seawater and excess water depth has no 
effect on burial depth. At depths from 0 to 20 meters, the variance in burial depth 
depends on both altitude and water depth since the vertical velocity of the mine as it 
enters the water becomes pertinent (Fig. 8). 

Water temperature was found to have no effect on the model’s outcome. 
Although temperature variance does alter the density of water up to 3% and also affects 
the viscosity (Stanley, 1969), this effect is not significant enough to alter the burial depth 
value calculated by the model. 

All cases discussed thus far assumed an initial angle of 90° with respect to the 
horizontal! and with a rotation rate of zero. This produces a situation where the mine is 
heading directly downward throughout the entire simulation, resulting in the maximum 
burial depth. When this initial attitude is varied, the burial depth is affected greatly as 


outlined in Table 3. 






a Altitude = 1.5 meters | Altitude = 150 meters 
Fall Angle = 0 
Fall Angle = 90° 


Table 3. Maximum and minimum burial values for a mine released from 1.5 or 150 
meters. An initial fall angle and subsequent sediment impact angle of 0° indicates a 
perpendicular orientation and maximum burial depth. Fall angle of 90° indicates the 
mine is parallel to the sediment and yields a minimum burial depth. 
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B. SENSITIVITY TO SEDIMENT CHARACTERISTICS 


As expected, sediment parameters are the most critical element in determining 
how deep the mine was buried when it came to rest. Sensitivity to the alteration of 
sediment density and shear strength was tested two ways. First, six sediment profiles 
were entered into the model and the resulting burial depth was examined (Fig. 9). These 
included three profiles from Sydney Harbor (Mulhearn, 1993) and three profiles available 
for selection in the IB model. The profiles included in the model are called simply 
“softsed”, “medsed”, and “hardsed” and do not clearly correspond to specific sediment 
types. Second, simplified cases of a single layer of sediment were used with constant 


density, varying shear strength and constant shear strength with varying density. 


1. Sensitivity to Shear Strength and Density 
Figure 10 illustrates the sensitivity of burial depth on density and shear strength. 


Here, a simple profile of just one layer was used and density and shear strength were 
varied separately. All other parameters were kept unchanged as default values. Plot (a) is 
| the burial depth with shear strength held constant and varying density from 1000 to 2000 
kg/m>. Shear strength of 1 kPa indicates extremely soft sediment, and density has a 
noticeable effect on burial depth of 37%. At the more common shear strength values of 5 
to 15 kPa density has little effect, just 3.7%. Plot (b) illustrates the effect of varying 
shear strength while keeping density constant. Again, we see the greatest ain of 
density value on the model output at low shear strength values. As shear strength 
increases, so does the influence of varying density. 

Tables 5 and 6 illustrate the effect on burial depth of simplifying the sediment 
profiles in two ways. Several methods for simplifying the sediment profile requirements 


were investigated, using the full profile case as a control. First, the density and shear 
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strength were held constant to 5 meters. The relative difference is under 26% for all 
profiles using this simplification. Next, a process of manufacturing sediment profiles 
using the density values measured at the water-sediment interface was derived and 
wie to the model. The profiles were assumed to consist of a constant value layer at 
the surface to a depth hj, a sharp gradient to ho, and then constant to a depth of 5 meters. 
The profiles were first applied to density only, holding shear strength constant, and then 
to both density and shear strength. Values for hj, ho, p(h2), and t(h2) were calculated by 
applying p, and t,. to polynomials derived from the data. The softsed and madsed 
profiles create the greatest differences from the control in all cases. 

Interestingly, creating a simulated density profile and keeping the shear strength 
value constant had no effect on the burial depth result when compared to keeping both 
values constant for five meters. This serves to underscore the fact that it is the shear 
strength parameter that has the primary influence on burial depth, not the more easily 


measured density parameter. 


2. Simulated Sediment Profiles 
Several attempts were made to manufacture shear strength profiles from density 


and shear strength values measured at the interface. This was explored in order to 
determine if a viable method of simplifying the data entry for the sediment phase could 
be devised. One attempt consisted of applying a fitted polynomial to measured density 
and shear strength values to create a synthetic profile from only interface values. Values 
for the sediment profiles used in the study, calculated with the following equations, are 


listed in table 5. 
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Based on the density profiles in figure 9, we empirically derived curve-fitting 


equations to represent the density profile, 


hl = -0.000061833*p, + 0.01609 

h2 = -0.0015*p, + 3.10 

S.(h2) = [(2.2048*10°)p,] — (0.109p*) + 201.1381p7[(1.6432*10°)p] + 
(5.0136*10’) 

p(h2) = Po + 250 


| softsed.sed_| medsed.sed_| hardsed.sed 
| bm) | 0.093 | 0.084 | 0.081 





Table 4. Values calculated using equations derived by fitting a polynomial to known 
density and shear strength profiles. The precision of this simplified method of depicting 
a sediment profile is demonstrated in table 5. 
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Mossman Rose B Off 
Bay *Y | Woolwich 
Measured Po, Suheld | 9 19) m 0.093 m 0.132 m 
constant to 5 meters 


Density and Shear 
strength profiles 0.101 m 0.059 m 0.105 m 
created using 
measured p, and S,, 


Measured po, Suheld | 9 693 m 0.463 m 0.094 m 
constant to 5 meters 
Density and Shear 
strength profiles 0.300 m 0.179 m 0.085 m 
created using 
measured p, and S,,, 


Table 5. Predicted burial depths using manufactured profiles based on measured 
values at the interface. 
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Hayter (1986) discussed an equation originally derived by Krone (1963) for 


deriving shear strength, S,, from density using empirically derived coefficients « and B: 
S.=ap* 
Values for a and B must be calculated for each separate sediment type, after which the 
shear strength can simply be calculated using the coefficients. F igure 11 illustrates the 
impact of varying a and B on the model output, given a constant density. The profile was 
assumed to consist of one layer of homogenous material. As expected, as a and B 
increase, shear strength also increases and burial depth decreases. Figure 12 is a series of 
contour plots with varied values as the axes. The contours represent predicted burial 
depth values. Known shear strength values are marked on the corresponding density plot. 
For all cases, there is a unique value of the coefficients that will produce a shear strength 
value given a specific density. Please sila that, while they are plotted here as one density 
value and one shear strength value per sediment type, a change in density would produce 


a corresponding change in shear strength that could be determined by use of the same two 


unique values of the « and B coefficients. 
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Figure 6. Derivative plot of all four output parameters from one IB 
model run. Notice the shape of each derivative plot is similar, 
confirming that one can be used to represent sensitivity levels of all 
output parameters. Burial depth was examined in the study for this 
paper, as it ts the primary result obtained by the model. 
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Figure 6. Continued from previous page. 
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Figure 8. Continued from previous page. 
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Figure 9. Sediment profiles of (a) density (kg/m*) and (b) shear strength 
(kPa) used in the sensitivity study. Data obtained from Mulhearn 


(1993). 
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Figure 10. Effect of (a)-density and (b) shear strength on burial depth. 
Density change only impacts the predicted burial depth in very soft 
sediments. As expected, shear strength has a dramatic impact on 
predicted burial depth. 
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Figure 11. Effect of a and f coefficients on predicted burial depth (m). 
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Figure 12. Contour plots of predicted burial depth (m) with 
respect to « and £ coefficients. Asterisks mark the values of a 
and § that correspond to the interface values of the three 
sediment profiles from Mulhearn (1993) used in this study and 
to a subset of the Voelkner (1973) data set. 
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VI. RECOMMENDATIONS 


As expected, there are several parameters that are both rarely known by the 
operator and of little import to the outcome of the model. In order to make the model 
easier to use with out sacrificing accuracy, these parameters should be simplified as much 
as possible. 

Water ‘alana was found to have no effect on burial depth, and should be 
eliminated from the list of variables. Altitude values should be simplified to represent a 
mine laying platform and the most likely height for release from that platform. For 
instance, it is more likely that an operator would know if the mines were laid by a ship or 
by a certain type of aircraft that the enemy has than at which altitude that aircraft was 
flying. Using intelligence and experience, one could form a rough database for selected 
platforms. The difference between an altitude of 1 meter and 300 meters is significant, 
but the difference between 1 meter and 5 meters is not. 

An equation for terminal velocity could be built into the data entry program that 
takes the weight of the chosen mine into account and asks the user if the water depth is 
less than that which would produce terminal velocity in the water phase. Assuming no 
rotation rate was chosen, if the depth was known to be greater than required for terminal 
velocity, a depth need not be entered. This simplification would also allow the user to 
have some confidence in the result as he moves about the area, regardless of water depth 
changes. 

The model could be revised to provide a range of values for burial depth, based on 
a initial attitude of 90° and 0°. In this way, the uncertainty of initial attitude and rotation 


rate would be eliminated and a more realistic range of values would be produced. This 
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may seem like a reduction in sophistication of the numerical model, but the reality is that 
the exact burial depth will never be known due to the unpredictability of the attitude of 
the mine as it encounters the sediment interface. This would also eliminate any effects 
due to currents or winds, since the primary effect of these influences would be on the 
attitude of the mine. 

Further investigation is warranted to simplify the sediment profile data entry. 
Assuming the values for density and shear strength are either known or can be measured 
at the interface, a set of equations should be derived and refined to create the remainder 
of the profile. If this were an option in the model, while still allowing the user to enter 
the entire profile 1f known, it would substantially increase the usefulness and precision of 
the model. The few cases discussed here and the equations derived from that limited data 
set are encouraging, and may indicate that such equations are possible and beneficial. 

A more precise method for simplifying the sediment data ingest would be to 
obtain samples from areas where mine clearing operations are likely to be accomplished. 
The relationship between density and shear strength for each sample would then be 
computed using the a and B coefficients discussed previously. The user could determine 
sediment type and density using available environmental atlases and the model would 


synthesize the profile. 
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